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Abstract 

— -^ , The sensitivity to the sterile neutrino mixing at very short basehne reactor neutrino experiments 

CNl ' is investigated. In the case of conventional (thermal neutron) reactors it is found that the sensitivity 

Q^' is lost for Am^ ^ 1 eV^ due to smearing of the reactor core size. On the other hand, in the case 

^ ■ of an experimental fast neutron reactor Joyo, because of its small size, sensitivity to sin 29i4 can 
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be as good as 0.03 for Am^ ~ several eV^ with the Bugey-like detector setup. 
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I. INTRODUCTION 

Since the LSND group announced the anomaly which suggests neutrino oscillations with 
mass squared difference of 0{1) eV^[lH3|, schemes with sterile neutrinos have attracted a 
lot of attention. This is because the standard three flavor scheme has only two independent 
mass squared differences, i.e., Am2i = Ami, ~ 8 x 10~^eV^ for the solar neutrino oscillation, 
and lArrig^l = Aml^^ ~ 2.4 x 10~^eV^ for the atmospheric neutrino oscillation, and it does 
not have room for the mass squared difference of (9(1) eV^. The extra state has to be sterile 
neutrino, which is singlet with respect to the gauge group of the Standard Model, because 
the number of weakly interacting light neutrinos has to be three from the LEP data, j^ 

The LSND anomaly has been tested by the MiniBooNE experiment. While the Mini- 
BooNE data on the neutrino mode [Si] disfavors the region suggested by LSND, their data 
on the anti-neutrino mode [6| seems to be consistent with the LSND data. 

Recently the flux of the reactor neutrino was recalculated in Ref. [7| and it was reported 
that the normalization is shifted by about +3% on average.^ Then a re-analysis of 19 reactor 
neutrino results at short baselines [9| in the light of new reactor neutrino spectra [7J may 
indicate neutrino oscillation at Am^ ^leV^. 

Reactor experiments with more than one detector have attracted much attention as a 



possibility to measure 6'i3 precisely |10l-ll3l|. and three experiments |14i4l6l| are now either 
running or expected to start soon. In the three flavor case with |Am|]^| = 2.4 x 10~^eV^, 
it was shown assuming infinite statistics that the optimized baseline lengths Lp and Lm 



for the far and near detectors are Lp ~1.8km and Ljy ~Okm in the rate analysis [17|, 118 
while they are Lp ^10. 6km and Ljv ~8.4km in the spectrum analysis [l9j. Unfortunately, 
in order to justify the assumption on negligible statistical errors for L ~10km, one would 
need unrealistically huge detectors, so one is forced to choose the baseline lengths which 
are optimized for the rate analysis for Am"^ = 2.4 x 10~^eV^. On the other hand, if we 
perform an oscillation experiment to probe Arn^ ~ (9(1) eV^, it becomes realistic to place 
the detectors at the baseline lengths which are optimized for the spectrum analysis (See 
Sect. 4 in the published version of Ref. |19|). 

In this paper we discuss the sensitivity of very short line reactor experiments to the sterile 
neutrino mixing for Arn^ ~ C^(l) eV^ in the so-called (3+l)-scheme.^ Proposals have been 
made to test the bound of the Bugey reactor experiment |24| on the sterile neutrino mixing 
angle using a reactor J25| ,^ an accelerator 26|, |27| , and a /3-source 28 



Throughout this paper we discuss the case with a single reactor and two detectors, which 
have an advantage to cancel the systematic errors that are correlated between detectors. The 
conditions of the detectors are assumed to be the same as those of the Bugey experiment. 

In Sect. 2 we briefly review the four neutrino schemes. In Sect. 3 we evaluate the sen- 
sitivity of very short line reactor experiments to the sterile neutrino mixing. In Sect. 4 we 
summarize our results. In Appendix |X] we give some details to get an analytical expression 



^ The increase in thfijlux claimed by Ref. [7| was confirnied by an independent calculation in Ref. [S]. 
^ According to Ref. [20|, the so-called (3+2)-schenie 21| gives a better fit to the global data than the (3+1 



scheme does. However, we are mainly interested here in the potential sensitivity to the sterile neutrino 
mixing for Arn^ ~ C'(l) eV^, and for simplicity we will discuss the (3+l)-scheme, leaving the analysis of 
the (3+2)-schenifi_for the future. 



•^ See, e.g., Refs. Il9| (the published version), 
reactor. 
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II. FOUR NEUTRINO SCHEMES 

Four-neutrino schemes consist of one extra sterile state in addition to the three weakly 
interacting ones. Depending on whether one or two mass eigenstate(s) are separated from the 
others by the largest mass-squared gap, the schemes are called (3+1)- and ('2+2)-schemes. 
The (2+2) schemes are excluded by the solar and atmospheric neutrino data (29|, so we 
will not discuss the (2+2) schemes in this paper. In the (3+1) schemes, on the other hand, 
the phenomenology of solar and atmospheric oscillations is approximately the same as that 
of the three flavor framework, so as far as the tension between the solar and atmospheric 
constraints are concerned, the (3+1) schemes do not have any problem. However, if we try to 
account for LSND and all other negative results of the short baseline experiments, then the 
(3+1) schemes have a problem. To explain the LSND data while satisfying the constraints 
from other disappearance experiments, the oscillation probabilities of the appearance and 



disappearance channels have to satisfy the following relation [30|, |31 



sin^ 26'LSND(Am^) < - sin^ 26'Bugcy(Am^) • sin^ 26'cDHsw(Am^) (1) 

where 6'LSND(Am^), ^cDHsw(Am^), 6'Bugey(Am^) are the value of the effective two-flavor mixing 
angle as a function of the mass squared difference Am^ in the allowed region for LSND 

fz7, 



— 7- z^e), the CDHSW experiment [32] [Vfj, -^ i^^), and the Bugey experiment |2J] (z/g — ^ K 



respectively. The reason that the (3+l)-scheme to explain LSND has been disfavored until 
Refs. [tI, IoJ appeared is because Eg. ([T]) is not satisfled for any value of Am^, if we adopt the 
allowed regions in Refs. |32j and |2J]. However, if the flux of the reactor neutrino is slightly 
larger than the one used in the Bugey analysis |2J|, the allowed region becomes slightly 
wider and we have more chance to satisfy Eq. ([1]). 

In t his p aper we will use the following parametrization for the mixing matrix, adopted 



in Ref. 33 : 



U = R^i{9'ii, 0) R2i{d2i-, 0) -R23(^23, ^3) -Rl4(6'l4, 0) -Rl3(6'l3, 82) Rui^u, Si) , 

where Rjk{Ojk, Si) are the complex rotation matrices in the jA;-plane deflned as: 

[Rjkidjk, Si)]pq = Spq + {C0s6jk - \){SjpS^q + 4p4g) + sin6'jfe(e"* '(5jp4g - e* '^jg^fcp)- 

With this parametrization, for the very short baseline reactor experiments, where the average 
neutrino energy E is approximately 4MeV and the baseline length is about 10m, we have 
I Am^^L/4_E| ^ 1 (j. A; = 1, 2, 3), so that the disappearance probability is given by 

P(^e ^ ^e) = 1 - sin^ 2^14 sin^ (^^) (2) 

to a good approximation. Eq. ([2]) is the formula which will be used in the oscillation analysis 
throughout this paper. 



III. SENSITIVITY TO sin^ 2014 BY A SPECTRAL ANALYSIS 

Throughout this paper we discuss the case with a single reactor and two detectors. The 
detectors are assumed to be of the Bugey type, i.e., hquid scintillation detector of volume 
600 liters with the detection efficiency which yields about 90,000 events at L=15m from a 
reactor of a power 2.8GW after running for 1800 hours. Also for simplicity, we assume in 
this paper that the near and far detectors are identical and have the same sizes of systematic 
errors. 

To evaluate the sensitivity to sin^ 26*14, we introduce the following y^ which was adopted 
in Ref . ET 
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Here, mf is the number of events to be measured at the near (A = iV) and far (A = F) 
for the i-th energy bin with the neutrino oscillation^, and if is the theoretical prediction 
without the oscillation, (a/^)^ is the uncorrelated error which consists of the statistical plus 
uncorrelated bin-to-bin systematic error: 
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where crj^ is the uncorrelated bin-to-bin systematic error. For simplicity we assume that 
sizes of the bin-to-bin uncorrelated systematic errors of the detectors and of the flux are 
independent of the energy. Also we take the choice of the bins in such a way that the number 
of events for each bin is equal: t^ = t"^ (i = 1, ■ ■ ■ ,n). a is a variable which corresponds 
to a common overall normalization error ctdb for the number of events, a^ (A = A^, F) 
is a variable which introduces the detector-speciflc uncertainties adB of the near and far 
detectors, a^ (i = 1, ■ ■ ■ , n) is a variable for an uncertainty aob of the theoretical prediction 
for each energy bin which is uncorrelated between different energy bins.^ a^j (A = A^, F) is a 
variable which introduces an energy calibration uncertainty cXcai and comes in the theoretical 
prediction in the form of (1 + (y.c^\)E instead of the observed energy E. Thus, the deviation 
vf (divided by the expected number of events) from the theoretical prediction tf due to this 



^ As we will see below, in the spectrum analysis, unlike in the rate one, it is not necessary for the near 
detector to have events without oscillation, since the dominant term in y^ looks at the difference between 
the maximum and minimum of the oscillation pattern. In fact, in Fig. [1] below there are regions in which 
i^N > i^F- Therefore, it may be misleading to call the detector at the shorter baseline length the near 
detector. Nevertheless, we use these words throughout this paper for simplicity. 

^ Here we follow the notation for the systematic errors in Ref. [19]. The first suffix of a stands for the 
property for the systematic error with respect to the detectors while the second is with respect to bins, 
and capital (small) letter stands for a correlated (uncorrelated) systematic error. The correspondence for 
the notation in Ref. |l2| is as follows: crdb = o-oxp, crdB = crfc, o-Db = cTshapo, ctdb = (^a- 
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In Eqs. (j4]) and ([5]), A'p is the number of target protons in the detector, T denotes 
the exposure time, La is the basehne for the detector A, F{E) is the flux of z/g, and 
cr{E) is the cross section of the inverse /3 decay Ue + p — ?■ n + e"*". ii^ is the en- 
ergy of the incident Ug and it is related to the positron energy Ee and the masses rrin, 
rrip of a neutron and a proton by i? = Eg + irtn — rrip = Eg + 1.3 MeV. E' is the 
measured positron energy and we will assume that the energy resolution is given by 
8%/\/E, i.e., R{Ee, E') = R{E — nin + nip, E') is a Gaussian function which describes 
the energy resolution and is given by R^E^, E') = {l/\/27ra)exp[—{Ee — E'Y/2a'^], where 
a = 0.08y/{Ee + me)/MeV = 0.08^y{E - 0.8/MeV)/MeV. Our strategy is to assume no 
oscillation in the theoretical prediction tf, to substitute the number of events with oscillation 
in mf-, and to see the sensitivity by looking at the value of x^- Since x^ is quadratical in the 
variables a, a^, oti, a^j, we can minimize with respect to these variables in Eq. (|3]) exactly.^ 
Since we assume no oscillation for the theoretical prediction tf, the quantity vf in Eq. (jl]) 
is independent of A{= N, F): 



.N 



Vi. 



(6) 



Here we will take the following assumptions for the systematic errors: 



C^db 


= 0.5%, 


0"dB 


= 0.5%, 


0"Db 


= 2%, 


C^DB 


= 3%, 


C"cal 


= 0.6%. 



(7) 

We will take the number n of bins n=32 and the energy interval 2.8MeV< Ey <7.8MeV. 
The measurement is supposed to continue for 1,800 hours, and the total numbers of events 
is expected to be approximately 90,000x(15m/LA)^ {A = N,F). 



In principle we could take a different convention, such as multiplying the measured numbers mf by the 
uncertainty (1 + a + a + a^), etc., and in fact it is done in some references. With such a convention, it 
becomes complicated to work with an analytical approach because x^ is not a Gaussian with respect to 
mf, which includes oscillation parameters, after the minimizations of a's. However, the difference between 
such a convention and ours affects only the higher orders in cr^'s. The conclusion on the sensitivity with 
such a convention should coincide with ours numerically. 



A. An analytical expression for x^ 

As was demonstrated in Refs. 18|, [l9|, it is instructive to have an analytic expression for 
X^, since it enables us to see which term becomes dominant in the spectrum analysis. Un- 
fortunately, it is difficult to evaluate x^ analytically with non-zero statistical errors and the 
energy calibration uncertainty (Tcai- In the limit of infinite statistic, an analytical expression 
for x^ can be obtained (cf. Eq. (10) in Ref. [19|])- On the other hand, in the presence of 
finite statistical errors (with equal numbers of events for each bin, tf = t^)^ we can obtain 
an analytical form in the limit of (Jcai = and the result is given by Eq. ( IA7|) in Appendix |X1 

In the present case, under the assumption that the total numbers of events A^'^ and A^^ 
are of order 10^, that the number n of bins is 32, and that the systematic errors are given 
by Eq. ([7]), we get 



sin^ 2^14 
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<Jn + {l/N^ + l/N^)/2 - {{1/N^ - l/N^y/i + (a^Jn)^}^^ ' 



(8) 
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where we have introduced the variables 
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and Uj are the orthonormal eigenvectors of some nxn unitary matrices (see Appendix |X]) . n 
is the number of bins and N"^ {A = N, F) stands for the total number of events {nxtf = N"^). 
C2 ^ cos 992, S2 = siny92, and ip2 are angles defined in Eq. I \A6\} to diagonalize a 2 x 2 matrix. 
Because CT'^y,/n ~ 1/{2N^), we have < 992 ~ tt/4. The numerators in all the four terms 
on RHS in Eq. ([H]) are divided by the number n of bins, because it is known in the case of 
infinite statistics [19] that {ui ■ (D^ ±D'^)}'^/n etc. are almost independent of n for n ^ 16. 
The first (last) two terms in Eq. (^ correspond to the rate (spectral) analysis, because 
ui oc (1, ■ ■ ■ , 1)^, while the other vectors Uj (j = 2, ■ ■ ■ ,n) are orthogonal to ui. The first 
(second) term looks at the sum (difference) of the total numbers of events at the far and 
near detectors. The dominant error in the ffist term is the systematic error a^^ which is 
correlated between the two detectors, while the dominant one in the second is the statistical 
error 1/N^ + 1/N^ plus the systematic one a^g which is uncorrelated between the two 
detectors. These first two terms are basically what was discussed in the rate analysis to 



study the sensitivity to sin^2^i3 in Ref. [ll]. The third (fourth) term examines the sum 
(difference) of the significance of the spectrum shapes of events at the far and near detectors. 
As statistic increase, the denominator of the fourth term becomes the smallest, and it is the 
fourth term which dominates ^ in the limit of infinite statistics. With the finite numbers 
of events, however, the last three terms do contribute to x^- 

In the limit of large mass squared difference Am|j^, the vector D^ becomes D^ — )■ 

(1/2)(1, ■ ■ ■ , 1), Ml ■ Z3^ ^ v^/2 and u^- ■ D^ ^ (j = 2, ■ ■ ■ , n). So the right hand side of 
Eq. dH]) becomes (4cr^g)~^, and the sensitivity in this limit is given by sin^ 26*14 ~ V^'^dbX'^ 
= V4 X 2.7x 0.03 ~ 0.1. Note that x^ = 2.7 gives a 90%CL bound on sin^ 2^14. 

B. Thermal neutron reactors 

A core of conventional (thermal neutron) reactors typically has a diameter of 4m, a height 
4m and its thermal power is typically a few GW. As we will see below, when we try to perform 
an oscillation experiment at very short baseline lengths, smearing of the reactor core size 
gives a nontrivial contribution to the results because the baseline lengths are comparable to 
the core size. 

We have computed x^ i^i Eq. ([3]) numerically in the case of a thermal neutron reactor. 
Fig. [1] shows the sensitivity to sin^ 26*14 as a function of the baseline lengths of the near and 
far detectors for Am4i = leV^, and 5eV^, respectively. We have assumed a thermal power of 
2.8GW, and a diameter of 4m and a height of 4m, and we have varied the baseline lengths in 



the range 14m^ L^ ^40m. Unlike the case of infinite statistics 19|, the statistical errors are 
important in the present setup of the detectors, and longer baseline lengths are disfavored. 

For Am4i = leV^ (Fig.[T]upper panel), the set {L^, Lp) — (17m, 23m) gives the optimum. 
In contrast to the rate analysis, in which the optimized baseline length of the near detector 
is LAr=Om to avoid oscillations, the spectrum analysis with {Ln,Lp) = (17m, 23m) looks 
at the difference between the maximum and minimum of the spectrum shape with neutrino 
oscillations at L^ and Lp mainly for the energy region E^, ~ 4MeV where the number of 
events are expected to be the largest. 

On the other hand, for Aml^ = 5eV^ (Fig- [H lower panel), the sensitivity to sin^ 2^14 is 
as low as 0.1 for most of the set of the baseline lengths {Li^,Lp). This implies that Am4i 
is so large that we have average over rapid oscillations for almost any baseline lengths and 
the only term which contributes in Eq. ([8]) is the first term, that corresponds to the rate 
analysis. 

Fig. |2] shows the sensitivity to sin^ 26*14 as a function of Ar?7,|i in the case of the baseline 
lengths {Lm,Lp) = (17m, 23m). For Am^i ^ 2eV^, it indicates that the sensitivity is no 
better than 0.1, which is basically the result of the rate analysis. The sensitivity in the case 
of a hypothetical point-like reactor, where all the conditions for the detectors are the same, 
is also given in Fig. [2] for comparison. Fig. [2] indicates that the sensitivity would be as good 
as several xlO~^ for a few eV^, if the core were point-like. So we can conclude that we have 
poor sensitivity for Am|i ^ 2eV^ because of the smearing effect of the finite core size of the 
reactor. We have computed the sensitivity also for different numbers of bins (n=16, 64), 
and verified explicitly that the result of Fig. [2] is almost independent of n for 16 < n < 64. 
This is because the statistical errors are dominant in x"^- 
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FIG. 1: The sensitivity to sin 26i4 as a function of the two baseUne lengths in the case of 
an ordinary thermal neutron reactor of a diameter 4m, a height 4m, a thermal power 2.8GW. 
Amli = leV^ (upper panel), 5eV^ (lower) is assumed. 
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FIG. 2: The sensitivity to sin 2^14. The red hne stands for the one for a thermal neutron reactor 
(of a diameter 4m, a height 4m, a thermal power 2.8GW) with two detectors at Lp = 23m and 
L]\j = 17m, while the blue line for the one for a hypothetical reactor with a point-like core. 



C. Experimental fast neutron reactors 

In the previous subsection, we have seen that the sensitivity to sin^ 26'i4 is lost because 
of the smearing effect of finite core size. Fast neutron reactors are known to have a high 
power density, and an oscillation experiments using them might be more advantageous than 
those with thermal neutron reactors. While experimental fast neutron reactors have a small 
power, it may allow experimentalists to put the detectors at a nearer location because of 
their experimental nature. 

Here we will discuss the concrete example of Joyo [34] with MK-III upgrade [35!], which 
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FIG. 3: The sensitivity to sin 2^14 as a function of the two basehne lengths in the case of Joyo 34] 
with MK-III upgrade [35] which is an experimental fast reactor (of a diameter 0.8m, a height 0.5m, 
a thermal power 0.14GW). A7n|^ = leV^ (upper panel), 5eV^ (lower) is assumed. 
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FIG. 4: The sensitivity to sin 20i4 of the Joyo reactor with two detectors. The Blue hne stands 
for the one in the case with two detectors at Lp = 4m and Ln = 3ni, while the red line for the 
one for in the case with two detectors at Lp = Sm and L^v = 4m. The sensitivity in the case of a 
conventional thermal neutron reactor (the green line) is also given for comparison. 

is an experimental fast breeder reactor.^ It has a relatively small size (a diameter of 0.8m, a 
height of 0.5m), and a relatively large thermal power 0.14GW. Its power density is approx- 
imately 500 kW/i, whereas a typical power density of thermal neutron reactors is approxi- 
mately 50 kW/i. 

Again we have computed x^ ^^ Eq. (^ numerically in the case of Joyo with MK-III 
upgrade. Fig. [3] shows the sensitivity to sin^ 29i4 as a function of the baseline lengths of 
the near and far detectors for Aml^ = leV^, and 5eV^, respectively. We have assumed the 



^ An experiment [365 "^^^ performed to detect neutrinos from a fast neutron reactor at Joyo, although they 
did not get sufficient statistical significance. 
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same conditions for the detectors (the size, the efficiency, the exposure time etc.) as those 
assumed in Fig. [H From Fig. [3] we see that the set {L^-, Lp) — (4m, 8m) gives the optimum 
for /S.m\^ = leV^ (Fig- E] upper panel). 

Fig. m shows the sensitivity to sin^ 26*14 as a function of Aml-^^ in the case of the sets 
of the basehne lengths (L^yLp) = (4m, 8m) and {Ln,Lp) = (3m, 4m). The set of the 
baseline lengths {L^, Lp) = (3m, 4m) has a peak for sin^ 2^14 ~ 4 eV^, while the other set 
{Li\i,Lp) = (4m, 8m) has sensitivity to sin^ 26*14 for a wider range of Aml^. The sensitivity 
in the case of a conventional thermal neutron reactor is also given in Fig. [2] for comparison, 
and the advantage of the case with Joyo is clear. Also in this case the statistical errors are 
dominant in x^, and the result of Fig. IHis almost independent of the number n of bins for 
16 < n < 64. 



IV. DISCUSSION AND CONCLUSION 

In the framework of the (3+l)-scheme, we studied the sensitivity to sin^ 26*14 of very short 
baseline reactor oscillation experiments by a spectrum analysis. The assumptions are that 
we have two detectors whose size and efficiency are exactly the same as those used at the 



Bugey experiment [24 . 



In the case of a conventional thermal neutron reactor, which has a core of a diameter 4m, 
a height 4m and a thermal power 2.8GW, by putting the detectors at L^ = 17m and Lp = 
23m, we obtain the sensitivity as good as several xlO~^ for Aml^ ^ leV^, but we lose the 
sensitivity above leV^ due to the smearing of the finite core size. 

In the case of an experimental fast neutron reactor with a core of a diameter of 0.8m, a 
height of 0.5m and a thermal power 0.14GW, on the other hand, we obtain the sensitivity 
as good as a several xlO^^ for leV^ ^ Amli < lOeV^ if we put the detectors at Ln = 4m 
and Lp = 8m. 

In both types of reactors with the Bugey-like detector setup, the statistical errors are 
dominant, and the sensitivity is almost independent of the number n of bins of the spectral 
analysis. The reason that the case of the experimental fast neutron reactor is competitive 
despite its small power is because the total numbers of events at L ~ several meters are 
comparable to those of the case with the thermal neutron reactor at L ~ a few x 10 meters. 

Since the best fit value obtained in Ref. J9| is Am4i ~ 2eV^, an experiment using an 
experimental fast neutron reactor offers a promising possibility. 

It should be emphasized that the flux uncertainty by 3% which was pointed out in Ref. [TJ 
does not cause a problem in our discussion. The reason is because the main contribution 
to x^ comes from the spectrum analysis, which looks at the difference of the maximum and 
the minimum of the energy spectrum and the effect of the overall normalization of the flux 
(the flrst term of the RHS in Eq. ([H])) gives little contribution to x^ iii the region of Am4i 
in which the sensitivity to sin^ 2^14 is better than 0.1. 

Although we may not be able to put the detectors at such a close location, we could 
increase the detector volume and/or the measurement period, so that the data size is the 
same as those described in this paper. While the experimental feasibility of the idea of 
putting detectors at a location very near to an experimental fast reactor is yet to be seen, 
it is worth investigating whether this method is technically possible. 
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Appendix A: Derivation of the covariance matrix 

By choosing the bins in such a way that tf = t^, Eq. ([3]) becomes 
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Let us redefine the variables 
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and let us introduce a vector notation 
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Following the discussions in the Appendix A in ^8| , the matrix element of the covariance 
matrix can be obtained as the expectation value of yfyf'- 

,AB 



{p)tr = (yfyf) 



^ ^f fda n [dy^ [da^ [ da^^fl f da, yfyf exp (-^ 

•^ A=N,F -J -J -J i=\-J V 

where the normalization M is defined in such a way that (1) = 1. From a straightforward 
calculation, we have 



Thus we have 



AB„.A„.A 2 



^i Vj ^cal- 



x' = {y^T^rnp''(lF 



in^20i4(D^^ D^T)p- 



sm 



D 



F / ' 



where the covariance matrix p is defined by 



P 



N 



M + l/t'ln 



(Al) 



with 



M = (^Db + ^db) In + (C^DB + (^Ib) Hn + (^l^lGr. 



N = alJn + alr^H, 



Dh^n -r UY^^±±n- 



(A2) 
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Here /„ is an n x n unit matrix, Hn and Gn are n x n matrices defined by 

1 ■■■ 1 



H„ 



(jn 



1 ■■• 1 

t;if2 vl 



\ 



VlVn V2Vn 



VlVn \ 
V2Vn 



(A3) 



(A4) 



where fj (j = 1, ■ ■ ■ ,n) is defined by Eqs. (jl]), and (|6]). Note that the covariance matrix 
does not include oscillation parameters but only errors. Diagonalization of the covariance 
matrix is useful to see which errors dominate x^- 

Unlike in Ref. [l9|, the presence of the statistical errors 1/t"^ and the matrix G„ makes it 
difficult to evaluate the eigenvalues analytically, so only in this appendix, we will take the 
limit al^i — t- which may be justified in the present setup with the total numbers of events 
~ 10,000 x(15m/LA)^. 

To diagonalize the matrix (jAip with a^^^^ = 0, we first note that the matrix Hn can be 
diagonahzed as 



H„ 



UnDnU-\ 



where 



is a diagonal matrix and 



Dn = diag(ra, 0, ■ ■ ■ , 0) 



Ml 



[Ul, 
1 



n 



U2 = ^(1,-1,0, 



.0)^ 



M3 



71^^'^' 



-2,0, 



.0)^ 



u„ = 



\/n(r 



l,-{n-l)y 



[n-1) 
is a n X n unitary matrix. Thus M and A^ can be diagonahzed as 

M = Un{ (4, + al^) In + (4b + ct^b) Dn} U-' 

= Undiag{pi,p2,--- ,P2)U-^, 

N = Un {alJn + al^Dn) U;^' 

= Undmg{qi,q2,--- ,(12) U;;^, 



14 



where 



Pi 


= ^Db + ^db + ^(4b + ^dB 


V2 


= ^Db + ^db 


Qi 


= (Job + nali^ 


q2 


— 2 
= ^Db- 



(A5) 



Now flAip can be cast into a block diagonal by 



a 



-1 







u 



-1 



Un 



u„ 



u 



M + l/t^In N 

N M + 1/t^In 

(^Db + ^db + V^"^) In + {ctIb + crl^) D„ al^In + al^D^ 

/fii ■■■ \ 

^2 ■•. : 







w 



V ■■■ B2J 

where U is a 2n x 2n unitary matrix whose matrix elements are given by 



U 



pq 



(I ■■■ 
■■■ 
1 






■•• 





■■■ 
^0 ■■• 


1 











•■■ 0\ 


1 





... 








... 





1 


... 

■•• 







1/ 



/ . {^p,23~l^q,3 + ^P,2j^q,n+j) ■ 

i=i 



pq 



Bi, B2 are 2x2 matrices defined as follows and are diagonalized as 



B2 



where 



qi 

qi pi + 1/t^ 

P2 + 1/t^ g2 

q2 P2 + l/t^ 



0-2 



(-) \(+h 



e*^^"Miag(Ap,ArOe"'^' 



;''^2"Miag(A^"\ A^+^) e-'^2"^ 



-i 
i 



is the Pauli matrix, (fi and ip2 sue introduced to diagonalize each matrix and are given by 

tan2(^. - ^g^ - 2Kb + ^^DB) 

tan^c^i - y^F_-^^^N- iHF_i/tN 

tan2<^2 = , ,^^' , ,, = , ^^"""V ,, , (A6) 
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.(±) 



and A^- (j = 1, 2) are the eigenvalues given by 



Ai*'- 



n 



n 



I ^2 



+ 



1 



1 1 
+ 



2 V A^^ N^ 



.(±) 




In the last expression, we have factored out n from the eigenvalues for later convenience, 
and we have introduced the total numbers of events A^"^ = n x t"^ (A = N,F). 
To evaluate x^, "we need to calculate the following factor: 



\ 

\ 
( 



■lip\(J2 












~«ip20"2 



-lVlO"2 








-l(,02O'2 



-l(/3lO-2 



-Ikp2cr2 



'%ip2CF2 



\ 



U2-D^ 

Un ■ Df 








-lip2CJ2 







"' Vf/„-. 






/ui-D''\ 








-tifi20-2 



u 



N 






J 



J 

\un-D^ J 

( ( ciui ■ D^ - siMi ■ D^ 
ciui ■ D^ + siui ■ D^ 
C2U2 ■ D^ - S2U2 ■ D^ 
C2U2 ■ D^ + S2U2 ■ 3^ 

C2Un ■ D^ - S2Un ' D^ 
\ V C2Un ■ D^ + S2Un " D^ 



where Cj = cos ipj, Sj = sinipj (j = 1,2). 
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Putting everything together, we obtain 

,2 



X 



sin^ 2^14 



(ciui ■ D^ + siui ■ D^\ In 



icxux ■ D^ - SiUi ■ D^\ /n 






fc2M,- ■ /3^ - S2U,- ■ Z)^) /n 

U (^L + O/n + (1/iV^ + l/iV^)/2 - {(l/ATi^ - l/N^f/^ + (^!,b/^)^}'/'' 

where both the numerators and denominators are divided by the number n of bins, because 
the numerators are almost independent of n after being divided by n in the hmit of infinite 
statistics 19 



In the present case, assuming that the total numbers of events A^^ and N^ are of order 
10^ and that the number n of bins is 32, our reference values (JTj) satisfy cr^g ~ 9 x 10"^ 3> 
l/N^ ~ 10-5 {A = A^, F), o-^B ~ 3 X 10-5 > al^/n ~ 1 x 10-^ alJn ~ 1 x 10"^ ~ l/iV^ ^ 
10-5 (^ = N,F). Hence we obtain Eq. ^. 
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